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Abstract 

We consider hard three-loop corrections to hyperfine splitting in muonium and positronium gen¬ 
erated by the diagrams with closed electron loops. There are six gauge-invariant sets of such 
diagrams that generate corrections of order ma^. The contributions of these diagrams are calcu¬ 
lated for an arbitrary electron-muon mass ratio without expansion in the small mass ratio. We 
obtain the formulae for contributions to hyperfine splitting that in the case of small mass ratio 
describe corrections for muonium and in the case of equal masses describe corrections for positro¬ 
nium. First few terms of the expansion of hard corrections in the small mass ratio were earlier 
calculated for muonium analytically. We check numerically that the new results coincide with the 
sum of the known terms of the expansion in the case of small mass ratio. In the case of equal masses 
we obtain hard nonlogarithmic corrections of order ma^ to hyperfine splitting in positronium. 
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I. INTRODUCTION 


For many years hyperfine splitting (HFS) in muonium and positronium remains an active 
field of experimental and theoretical research. Results of highly accurate HFS measurements 
can be compared with the theoretical predictions of quantum electrodynamics obtained from 
the first principles without any adjustable parameters. Both experiment and theory have 


achieved very hig’ 
interval 16-51 Hz 


1 accuracy. The experimental errors for HFS in muonium are now in the 
,l2|], and a new measurement with the goal to reduce the error to about 10 




Hz or to a few parts in 10® is now planned at J-PARC |3|, 1^. Current theoretical uncertainty 
of HFS in muonium is about 70-100 Hz, see, e.g., reviews in js-?]. Recent theoretical work 
on HFS in muonium concentrated on calculation of radiative-recoil corrections of order 
m/M)Ep that arise from the three-loop diagrams with closed electron and muon loops 



'm 

8- 
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12l |. The goal of this work is to reduce the theoretical error below 10 Hz. 


The hyperfine splitting in positronium is measured with the error bars at the level of 1-2 


MHz 
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16l |. There is a discrepancy about three standard deviations between the results 


of old and new experiments. New measurement of the positronium HFS splitting is now 
planned at J-PARC 17|. All theoretical contributions to HFS in positronium of order ma^ 
and logarithmic corrections of order ma^ are already known, see, e.g., reviews in 181-120 1. 


A new stage in the theory of positronium HFS was opened in 1^ where the one-photon 
annihilation contribution of order ma^ was calculated. This paper was soon followed by the 
works of Adkins and collaborators 20|, l2l|, who calculated contributions of the light-by-light 


scattering insertion in the scattering and annihilation channels. 

Hard nonlogarithmic contributions to HFS in positronium of order ma^ are similar to 
the radiative and radiative-recoil corrections to HFS in muonium of orders a^{Za)Ep and 



a 
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12| 


12l |. The corrections in muonium are power series in the electron-muon mass ratio 


with the coefficients enhanced by large logarithms of this mass ratio. The goal of the old work 
on muonium was to calculate the coefficients in this expansion, at least the factors before 
the logarithms, analytically. In the case of positronium the masses are equal and the hard 
corrections of order ma^ are pure numbers. We apply the approach developed for muonium 
to positronium. We consider an electromagnetically bound system of two particles with 
arbitrary masses M and m, and obtain general expressions for the hard corrections to HFS 
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of order ma^ without expansion in the mass ratio of the constituents. We check numerically 
that in the case of a small mass ratio these formulae reproduce with high accuracy the 
sum of all already known terms in the expansion in the small mass ratio for muonium. We 
use the general expressions for the case of equal masses and calculate all hard three-loop 
contributions to HFS in positronium of order ma^ that are due to the diagrams with closed 
electron loops. The results of these calculations were reported in the rapid communication 


22l |. Below we present the details of the calculations in the general case of arbitrary mass 


ratio and in the special case of equal masses, for positronium. 


II. CALCULATIONS 


We start with the infrared divergent contribution to HFS in muonium generated by the 
two-photon exchange diagrams in Fig. [1] calculated in the scattering approximation 


AE = 


Za 3mM 
- Ep- 


71 16 
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■^^L^ZkeM)Ll.,skel,a0{-q) 


TT 




nr'^q^ 


( 1 ) 


where 


jaP 

'^e.skel 


(?) = 


2g2 


-yqY = 2L2ez7^?7 




( 2 ) 


is the forward electron Compton scattering amplitude in the tree approximation (the skeleton 
electron-line factor), and is a similar amplitude for the muon. The Fermi energy is 

dehned as Ep = {8/3){ZaYml/{mM), where = mM/{m + M) is the reduced mass. In 
the case of equal masses, M = m, the Fermi energy Ep turns into the leading nonannihilation 
contribution to HFS in positronium Ep^ = ma'^/3. The external electron and muon lines in 
the diagrams in Fig. [T]are on the mass shell and carry zero spatial momenta. In the second 
line in Eq. ([1]) we calculated projection of the matrix elements on HFS. 
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FIG. 1. Diagrams with two-photon exchanges 
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After the Wick rotation and transition to four-dimensional spherical coordinates (go = 
qcos9, |g| = gsin6') we obtain 


ae = 

71 71 

Za AmM 
=- Ep - 

TT 71 


desiv?e / dg2(2 + cos2 0)4"2ei4tei 

Jo 

2 -I- cos^ 6 


dOsin^e 


dq^ 


'0 


(g2 + Am‘^ cos^ 9){q'^ + cos^ 6) 


(3) 


_ Za mM 

71 ^ 


dq'^ffiiq), 


where at the last step we rescaled the integration momentum g qm. The dimensionless 
weight function /^(g) in terms of an auxiliary function 


has the form 


fij) 


1 \/¥Ta 2^^^ 

4 4g g^ 


(4) 


fM) = Hq) -4/iV(2/ig), 

where /i = m/(2M). 

In the case of positronium M ^ m and the weight function simplihes 


(5) 


mM lQ + 2q^ + q'^ - q^^Jq^ + 4 
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Respectively, the skeleton integral in Eq. ([3]) in the case of positronium turns into 
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The integrals in Eq. (|3]) and Eq. ([7]) are sums of an infrared linearly divergent integral 
and a hnite one. In a more accurate approximation (with the off mass shell external fermion 
lines) the linear divergence is cutoff at the characteristic atomic scale ~ ma and turns into a 
contribution of a lower order in a. The remaining hnite part of the integral originates at hard 
integration momenta ~ m (or in the interval from m to M in the case of unequal masses) 
and generates a contribution of order aEp. In the case of unequal masses, for muonium, the 
linearly infrared divergent contribution turns into the leading nonrecoil Fermi contribution 
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Ep to HF^ while the hnite part generates the leading recoil correction of order a{m/M)Ep, 
see, e.g., p, l6|. Let us emphasize that due to the linear (as opposed to logarithmic) nature 
of the apparent divergence it leaves no hnite remnant of order aEp and should be simply 
thrown away. No need in matching of high and low integration momenta arises. 

Six gauge-invariant sets of diagrams in Figs. [2]-[Hand in Figs.El-IHlgenerate hard radiative 


corrections of order that are due to the graphs with closed electron loopi^ . All these 
diagrams can be interpreted as the results of radiative insertions in the skeleton diagrams 
with two-photon exchanges in Fig. [H It is well known that insertion of radiative corrections 


the low integtation .nomentum tegton, eee. e.g., Hence, all diagtants in 

Figs. [2]-[Hand in Figs. 0-El are infrared convergento. Moreover, the characteristic integration 
momenta in these diagrams are hard (of order ~ m or in the interval from m to M in the 
case of unequal masses) and are much larger than the atomic momenta of order ~ ma, what 
justihes validity of the scattering approximation for their calculation. This is exactly the 
approximation we used above in calculation of the contribution of the skeleton diagrams in 
Fig. [H and all corrections calculated below are obtained by some modihcations of the basic 
integrals in Eq. ([3]) and Eq. (171) . 


A. Analytic Results for One- and Two-Loop Polarization Insertions in the Ex¬ 
changed Photons 
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FIG. 2. Diagrams with two one-loop polarization insertions 

Consider hrst the diagrams in Fig. |2]with two one-loop polarization loops. Insertion of a 
polarization operator in a photon line with momentum q (all momenta below are measured 
in units of the electron mass) reduces to the replacement in the photon propagator 

4 ^ ( 8 ) 

TT 

^ All gauge-invariant sets of diagrams include the graphs with the crossed exchanged photons that we do 
not show explicitly. 

^ Linearly infrared divergent contributions due to the anomalous magnetic moment should be subtracted 


from radiative corrections in Figs. I6] and Fig.ISl sge more on this below. 















where (a/7r)Ji(g) is the well known representation of the one-loop vacuum polarization 241 


-h{q) = - / 
TT TT Jo 




1-^.2 ^2 , 4 


(9) 


r+ 13 

We see that a photon line that carries a polarization loop has a natural interpretation as a 
propagator of a massive photon with mass squared = 4/(1 — n^). According to Eq. ([9]) 
this propagator should be integrated over v with the weight (a/7r)n^(l — n^/3)/(l — n^). 

The contribution of the diagrams in Fig. [2] is obtained by insertion of the one-loop photon 
polarization squared (a/7r)^g^/f (g) in the integrands in Eq. (|3]) and Eq. ([7]). Due to nonsin¬ 
gular behavior of the polarization operator at —?■ 0 we obtain convergent integrals where 
the effective integration momenta are hard, of order ~ m (or in the interval from ~ m to M 
in the case of unequal masses). Then in the general case of unequal masses the contribution 
to HFS has the form 


AE = 3 


a'^{Za] 


E, 


mM 


d<l fMl A (9). 
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TT"' 

where the factor 3 before the integral has the combinatorial origin. We checked numerically 
that in the small mass ratio limit this integral reproduces the sum of all known analytically 
terms 2^, [2^ of the expansion of this contribution in the small mass ratio. 


In the case of positronium the integral in Eq. ffTOj) reduces to (compare Eq. ([7j)) 


AEi = 3^E^/ / dq^U{q)q^P,{q), (11) 

Jo 

and after computation we obtain the contribution to HFS of the diagrams with two one-loop 
polarization insertions in Fig. |2] 


AEi 
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FIG. 3. Diagrams with two-loop polarization insertions 
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The contribution of the two-loop vacuum polarization in Fi^ 


insertion of the two-loop photon polarization /Ti‘^)q^l 2 {q) 
Eq. ([2D and Eq. (|7D 


24 


3] ca n be obtained by the 


27( 1 in the integrands in 
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+ -n(5-3n2) I, 


(13) 


where the dilogarithm is defined as Li 2 (^) = — dtln(l — zt)/t. 

In the case of unequal masses the integral for the contribution to HFS of the diagrams 
with the two-loop polarization in Fig. [3] loop has the form 


. ^ a^(Za) ^ mM 

AE = 2 — 


7 |-o ]\/[‘2 _ 17^2 


dq^ff,{q)q‘^l2{q), 


(14) 


where the factor 2 before the integral is due to combinatorics. Again, due to nonsingular 
behavior of the two-loop polarization at small —?• 0 the integral in Eq. flTTl) is convergent, 
and typical integration momenta are hard, in the interval from m to M. We checked nu¬ 
merically that in the small mass ratio case the integral in Eq. (ITTp coincides with the sum 


25 


26| of the expansion of this contribution to HFS in the small mass 


of the known terms 
ratio. 

In the case of equal masses, for positronium, the contribution to HFS of the diagrams in 
Fig. [3] reduces to the integral 


a 


AE2 = 2^Ep / dq^U{q)q^h{q). 
TT'* Jo 


This integral admits an analytic calculation, and we obtain 


(15) 


AE 2 = 
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<3 +- n2 +-+ 

^ 15 675 360 


—E])" = 5.209 219 6U—Ep. (16) 
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B. One-Loop Electron Factor and One-Loop Polarization Insertion in the Ex¬ 


changed Photon 


.Z22. 






FIG. 4. Diagrams with one-loop polarization and radiative photon insertions 


The diagrams in Fig. 0] are obtained from the skeleton diagrams in Fig. [T] by one-loop 
radiative insertions in one of the exchanged photons and one of the fermion linej^. To 
describe these radiative insertions it is convenient to introduce the one-loop electron factor 
that is defined as a gauge invariant sum of the diagrams in Fig. [5] where the external electron 
lines are on-shell and carry zero spatial momenta (plus the diagrams with the exchanged 
external photon vertices). Physically the electron factor is a sum of one-loop corrections to 
the spin-dependent amplitude of the virtual forward Compton scattering. 

The gauge invariant electron factor can be written as a sum of two gauge invariant 
terms , where the tenn is the contribution of the anomalous magnetic 




moment (for more details see, e.g., 28|, l29|). The multiloop electron factors also can be 


written as sums of two gauge invariant terms. Representation of the electron factor in the 
form of a sum of two gauge invariant terms is convenient for calculations because these 
terms have different behavior at low virtual photon momenta. According to the generalized 
low-energy theorem (see, e.g., js, 6|) all terms linear in the small photon momentum q 
are due to the term while the term decreases at least as at small q^. This 
different low-energy behavior determines the structure of the integrals for the contributions 
to hyperfine splitting. In the case of the diagrams in Fig. |6]and in Fig. [Hlthe contributions 
to HFS generated by the term are of lower order in a than the apparent order of a 
diagram. Technically presence of the previous order contribution reveals itself as a linear 
infrared divergence of an integral calculated in the scattering approximation. 

In the diagrams in Fig. IHthe skeleton fermion line in Fig. [T]is effectively replaced by the 
one-loop fermion factor in Fig. |5l what can be described by the substitution 


Multiplicity factors in these diagrams correspond to the case of positronium, not muonium. 
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FIG. 5. One-loop fermion factor 


LZkeii^) ^ = 2 ^ { ql) + go 
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(17) 


where Li(n) are scalar form factors. The scalar form factors L^n) have the form 


Li — Li + La, 


Lit — Lit — L. 


(18) 


where the scalar form factors L^nj and La correspond to L^j, and L^“\ respectively. The 
factor 2 before the braces arises becanse we normalize the scalar form factors like the skeleton 
one in Eq. (|2]), and the factor Q;/(47r) is dne to the one-loop integration in the fermion factor. 

The one-loop electron factor L^j, with the snbtracted contribntion of the anomalons 
magnetic moment enters calcnlations of the two-loop radiative-recoil corrections to HFS in 
mnoninm, and we had derived an explicit integral representations for the respective scalar 
form factors L^n) long time ago 3nl-l32| . After the Wick rotation, rescaling of the integration 
momentnm q —)■ qm, and transition to the fonr-dimensional spherical coordinates the form 
factors can be written as 


r / 2 2a\ A r A ] (g^ + «^) [(g^ + - 126^^2 COS^6'] . 2 • , 4\ 
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where = x'^/y{l — y), b = {1 — x)/{1 — y), and the coefficient fnnctions c* are collected in 
Table E 
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TABLE I. Coefficients in the Electron-Line Factor 
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C3 

C4 
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C6 
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1 - 6x - 2x2 - f f26 - ^ - 37x - 2x2 -h 12xy -h 161nx 


The scalar form factor La is proportional to the respective skeleton form factor Lg^ei in 
Eq. (j2]). After the Wick rotation and in terms of the dimensionless integration momentum 
q it has the form 


La — 2 I A 2 a ~ ‘^Lskel 

-\- A cos^ o 

In the case of unequal masses an analytic expression for the contribution to HFS of the 
diagrams in Fig. IHis obtained by modihcation of the skeleton integral in Eq. ([3]). First, we 
replace the skeleton factor in the integrand 



(2 -1- cos^ 9)L 


(e) 

skel 


-)■ 
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dvr 


(2 + cos^ 6')Li — 3 cos^ dL^ 


( 21 ) 


The factor a/{A'K) comes from the substitution in Eq. (1T7|1 . and the term with arises 
because the one-loop electron factor in Eq. flT7|) contains an additional spinor structure in 
comparison with the skeleton one in Eq. (|2]). 

Second, we need to account for the polarization loops in Fig. 0] and insert the term 
2q^Ii{q) in the integrand in Eq. ([7]). The factor 2 is due to two ways to insert the polarization 
operator in one of the exchanged photons. The polarization operator q^Ii{q) decreases like 
at small g, and we obtain an infrared convergent integral with hard characteristic integration 
momenta of order m (or in the interval from m to M in the case of unequal masses). Due 
to suppression of the small integration momenta the anomalous magnetic moment in the 
diagrams in Fig. 0] gives contribution on par with the other terms in the one-loop electron 
factor 
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Then the contribution to HFS of the diagrams in Fig. 0] in the case of unequal masses 
has the form 


, ^ a^iZa) ^ M2 
AE = ^ ^ Ep 
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The leading terms of the expansion of the contribution to HFS in Eq. (|2^ in the small 


mass ratio are already known for some time 
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3 2 9 2 45 9 8 675 ^ 

m/5,nM 10, M \ 
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(23) 


We have checked numerically that the integral in Eq. (12^ coincides with this analytical 
result in the case of small mass ratio. 

In the case of equal masses there is an extra factor 2 before the diagrams in Fig. 01 This 
factor arises because now there are two ways to insert the fermion factor in one of the lepton 
lines. Hence, the respective contribution to HFS is described by the doubled integral in 
Eq. fl2^ at M = m. Then we obtain the contribution of the diagrams in Fig. 0]to HFS in 
positronium in the form 
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= -1.287 09 (1) ^E^\ 
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C. One-Loop Polarization Insertion in the Electron Factor 



FIG. 6. Diagrams with one-loop polarization insertions in radiative photons 
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Consider now the diagrams in Fig. [6] with the one-loop polarization insertions in the 
radiative photoru. Effectively these diagrams contain a massive radiative photon, see Eq. ([8]) 
and Eq. ([2]). In principle, the respective electron factor can be obtained from the one-loop 
electron factor in Eq. ffT7|l by restoring the radiative photon mass sqnared = 4m^/(l — 
n^), followed by integration over v with the weight (a/7r)n^(l — n^/3)/(l — n^). However, 
the relatively compact expression in Eq. (1T7|) is a resnlt of nnmerous cancellations in the 
integrand between the contribntions from different diagrams in Fig. |5l and technically it 
is mnch easier to start calcnlation of the two-loop electron factor in Fig. |6] from scratch. 
We consider this electron factor as a snm of the contribntions corresponding to the separate 
diagrams in Fig. Elwith the self-energy, vertex and spanning photon insertions in the electron 
line. Each of these terms is calculated as a one-loop diagrams with a massive photon and 
then integrated over the auxiliary parameter v as we just explained. 

The only subtlety in further calculations is connected with the diagrams with the vertex 
correction in Fig. [6l All entries in the two-loop fermion factor except the two-loop anomalous 
magnetic moment carry at least one extra power of at —)■ 0 in comparison with the 
skeleton electron factor. One can separate the contribution to the two-loop anomalous 
magnetic moment from the two-loop vertex in the second diagram in Fig. [6] in a gauge 
invariant way, like we separated the one-loop anomalous magnetic moment from the one- 
loop electron factor above. The two-loop anomalous magnetic moment term in the 
second diagram in Fig. [6] generates a linearly infrared divergent contribution to HFS. This 
linear infrared divergence that is cutoff at the characteristic atomic scale ~ ma indicates 
that the anomalous magnetic moment generates a contribution to HFS of the previous order 
in a. This correction of order ma^ is already accounted for in earlier calculations and 
we should simply delete the apparently divergent term that generates it. To get rid of 
the spurious divergence we subtract the gauge invariant term with the two-loop anomalous 
magnetic moment from the two-loop electron vertex in Fig. |6l The subtracted two-loop 
electron factor in Fig. [6] can be written in terms of scalar two-loop form factors Li „ 


r (2) _ j (2,S) o r (2,A) . y (2,=.) 

-^1,11 - -^1,11 “T II Ui II 


(25) 


similar to the one-loop form factors Li n in Eq. (ITSll . Unlike the one-loop form factors Li n in 


Multiplicity factors in these diagrams correspond to the case of positronium, not muonium. 
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Eq. (fT7|) these two-loop form factors do not include contributions of the anomalous magnetic 
moment. 

We have derived explicit expressions for the two-loop form factor in Eq. 025 p in 
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where we calculated nonrecoil and radiative-recoil correction to HFS in muonium due to the 
diagrams in Fig. O These expressions are rather cumbersome and we will not reproduce 
them here. The contribution to HFS of the diagrams in Fig. [ 6 ] in the case of unequal masses 
has the form 


AE = ^ Ep—— [ dq^ [ d9 siri^ 9 {2 + cos^ 9)l['^^ — 3 cos^ 9L[i^ . (26) 
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m TT 


A few leading terms of the expansion of this contribution to HFS in the small mass ratio 
were calculated earlier 
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TT"’ 1 ^ M 

We have checked numerically that the expression in Eq. fl26P derived for arbitrary masses 
reproduces the expansion above in the case of small mass ratio. 

In the case of equal masses the contribution of the diagrams in Fig. [ 6 ] to HFS reduces to 


AF4 = f 

TT'* TT Jo 


dq^ / d9 sin^ (^L^skei + cos^ 9)Lf^ — 3 cos^ 9L 


( 2 )' 
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where an extra factor 2 before the integral (in comparison with Eq. fl26l) ) arises because we 
can insert the two-loop electron factor in either of the fermion lines. 

After calculations we obtain 
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AF 4 = -3.154 41 (1) —E 


7lPs 




(29) 


D. Light-by-Light Scattering Contribution 

Due to gauge invariance the light-by-light scattering block fast decreases with the mo¬ 
menta of the external (virtual) photons. Therefore, effectively all integration momenta in 
the diagrams in Fig. Hare hard, of order of the electron mass (or in the interval from m to 
M in the case of unequal masses). 
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FIG. 7. Diagrams with light-by-light scattering insertions 


The contribution of the light-by-light scattering block to HFS in the general case of 
unequal masses has the form 12|, [29 1 


TT^ 327r Jq 


de sin^ 9 


T{q^^ cos^ 9) 
m?q^ + AM^ cos^ 9 


(30) 


The explicit integral representation for the function T(g^, cos^ 9) can be found in 12|. 

The first terms of the expansion of this contribution to HFS in the small mass ratio were 
calculated during the years Q Q, y, Q 
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We have checked numerically that the general expression in Eq. (I30|) coincides with the sum 
in Eq. fl3T]) in the case of small mass ratio. 

In the case of equal masses the integral in Eq. (I30|) reduces tcQ 


r 2COS^ 0) 
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^ 71^ ^ 327r Jq Jo ^ + A cos^ 9 


j d9 sin^ 9 
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Calculation of this contribution to HFS in positronium proceeds exactly like ca! 
the respective nonlogarithmic radiative-recoil correction to HFS in muonium in 
obtain 


(32) 

culation of 
and we 
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AE 5 = -0.706 27 (5) -^Ep, 
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what coincides with the result hrst obtained in 


20 |. 


There is a misprint in the respective expression in Eq.(14) in 


22 |. 


( 33 ) 
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E. Two-One Loop Electron Factors 



FIG. 8. Diagrams with one-loop radiative photon insertions in both fermion lines 

The diagrams in Fig. [S] contain the one-loop fermion factors from Eq. fll7p in both fermion 
lines. Naively, the contribution of these diagrams to HFS can be obtained from the skeleton 
integral in Eq. (|3]) by the replacement 


- (e) r (/i) 


/ a 


^skel^skeA‘^ + COS 6 ) -)■ — j 


(2 + cos^ - 3 cos^ e ( + Wl\ 


he) r (m) 


'(e) T ifJ.) 


(e) r (m) 


-f COS 6{1 + 2 cos 6 )Lu’L 


(34) 


where the terms on the right hand side arise after calculation of the projection of the product 
of two electron factors (see Eq. (flTj) ) on the HFS structure. 

The scalar form factors include terms with the scalar form factors arising 

due anomalous magnetic moments, see Eq. (ITSD . Therefore, each product of the scalar func¬ 
tions in the square brackets on the right hand side of Eq. fIMl) contains the term 
As we already mentioned the form factors decrease at least as at —?• 0 relative 

to the skeleton form factors, while the form factors behave exactly like the skeleton 

form factors, see Eq. fl 2 Up . Hence, each integral of is a sum of a linearly infrared 

divergent and finite contributions, compare with the skeleton integral in Eq. ([3]). In a 
more accurate calculation the linearly infrared divergent contribution would be cutoff at 
the atomic scale ~ ma and would generate a correction of lower order in a. It should be 
simply subtracted, while we need to preserve the hnite part of the integral that generates 
correction of order ma^. In the general case of different masses (for example, for muo- 
nium) the finite part is a recoil contribution and it was calculated in 2^. It is equal to 
(9/16)(mM)/(M^ — m^)ln(M^/m^) up to a normalization factor. Hence, subtraction of 
the linearly infrared divergent contribution due to terms with in Eq. (IMP reduces 

to addition of (9/16)(mM)/(M^ — m^)ln(M^/m^) to the respective contribution to HFS 
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instead of all terms on the right hand side in Eq. flM|l that proportional to After 

this replacement of the infrared divergent part we obtain a convergent integral with hard 
characteristic integration momenta in the interval from m to M. The contribution to HFS 
of the diagrams in Fig. |H]in the case of unequal masses has the form 


AE = J £ J_ 

TT^ \ m 4:71 


/ dq^ / dO sin^ 9d9 
'0 Jo 


(2 + cos^ 9) 


X ( + if 


- 3 cos" 9 ( - L'f'Lr + L 


(e) r (m) I r (e) r (m) 


+ cos" 0(1 + 2 cos" 9) 


+ 


9 mM 


In 


M" 


16 — m? m? 


(35) 


The hrst terms of the expansion of the contribution to HFS of the diagrams in Fig. [8] in the 
small mass ratio are already known for some time 28|, l37| 


AF 


a{Z‘^a){Za) ^ 

-5- Ef 

m f 9, M 

— — In- 

M \ 8 m 


TT 

Y 


In 2 


13 

T 


15,. . IStt" 

—(3) H--— In 2 + 


377r" 

24 


147\ 9 Mm m2 

32 / 16 M^ — m2 m2 


(36) 


We have checked numerically that in the case of small mass ratio the expression in Eq. fl3^ 
coincides with sum above. 

In the case of equal masses the integral in Eq. (l35ll simplihes 


AEq = %E 


Ps 


'•OO /•TT 

,2 


TT' 


^ 1 dvr 


dq / d9 sin 9d9 


{2 + cos^9){Y,+2LaL,) 


— 6 cos2 0(LiL„ + LaLii — LaLJ) + cos2 0(1 + 2 cos2 0)(L"j — 2LaLii) 


(37) 


+ 


16 


After numerical calculations we obtain contribution of the diagrams in Fig. [S] to HFS in 
positronium 


AEg = -4.739 55 (40) 


-^Ep 


Ps 


(38) 
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III. SUMMARY OF RESULTS 


We have derived explicit expressions for hard three-loop contributions to hyperfine split¬ 
ting generated by the six gauge-invariant sets of diagrams with closed electron loops in 
Figs. [2]- [Hand in Figs. [6]- 0 In the case of unequal lepton masses we confirm numerically 
the already known results for muonium obtained earlier in the form of an expansion in the 
small mass ratio. We have calculated the contributions of these diagrams to HFS in the 
case of equal masses, for positronium. Collecting the results in Eq. ([12]), Eq. (fT6|) . Eq. fl2^ . 
Eq. ((29|), Eq. fl33|) . and Eq. fl38|) . we obtain the total hard contribution to HFS in positron¬ 
ium of order generated by all diagrams with closed electron loops in Figs. [2]- Hjand in 
Figs. El - 0 


a\ 3 


AE = -3.875 0 (4) ( - ) = -1.291 7 (1)11^ = -5.672 kHz. (39) 


ma 


TT 


TT"’ 


Taking into account all other recent theoretical results 19l-l2l[| we obtain the theoretical 
prediction for HFS in positronium 


AEtheor = 203 391.89 (25) MHz. 


( 40 ) 


The latest experimental result is 



AEe,p = 203 394.2 (1.6),*^* (1.3),^,, MHz. 


(41) 


There are no contradictions between theory and experiment at the present level of accuracy, 
but further reduction of both the experimental and theoretical uncertainties is warranted. 
Calculation of the remaining ultrasoft and hard nonlogarithmic contributions of order ma^ is 
the next task for the theory. We hope to report the results for the remaining hard corrections 
in the near future. 
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